Comparing the temperatures of galaxy clusters from hydro-N-body simulations to 

Chandra and XMNl-Newton observations 



P. Mazzotta^, E. Rasia^, L. Moscardini^, G. Tormen^ 
^Dip. di Fisica, Universitd di "Tor Vergata", via della Ricerca Scientifica 1, 1-00133 Roma, Raly 
Dipartimento di Astronomia, Universitd di Padova, vicolo deWOsservatorio 2, 1-35122 Padova, Raly 



Theoretical studies of the physical processes in clusters of galaxies are mainly based on the 
results of numerical simulations, which in turn are often directly compared to X-ray obser- 
vations. Although trivial in principle, these comparisons are not always simple. We show 
that the projected spectroscopic temperature of clusters obtained from X-ray observations is 
always lower than the emission-weighed temperature. This bias is related to the fact that the 
emission- weighted temperature does not reflect the actual spectral properties of the observed 
source. This has implications for the study of thermal structures in clusters, especially when 
strong temperature gradients, like shock fronts, are present. In real observations shock fronts 
appear much weaker than what is predicted by emission-weighted temperature maps. We 
propose a new formula, the spectroscopic-like temperature function that better approximates 
the spectroscopic temperature, making simulations more directly comparable to observations. 



1 Introduction 

Recent observational data with high spatial and spectral resolution suggest that clusters are far 
from isothermal, and show instead a number of peculiar thermal features, like cold fronts 
cavities , blobs and filaments . 

A theoretical interpretation of these observations clearly requires state-of-the-art hydro-N- 
body simulations, which can be used to extract realistic temperature maps and/or profiles. The 
comparison between real and simulated data, however, is complicated by different problems, 
produced both by projection effects and by instrumental artifacts. A further complication can 
arise from a possible mismatch between the spectroscopic temperature Tgpec estimated from X- 
ray observations and the temperatures usually defined in numerical results. In fact, while Tgpec is 
a mean projected temperature obtained by fitting a single or multi-temperature thermal model 
to the observed photon spectrum, theoretical models fully exploit the three-dimensional thermal 



Figure 1: Comparison of the map of the emission- weighted temperature Tew (left-panel) and the spectroscopic- 
like temperature Tsi (right-panel) for the simulated cluster of galaxy. Both maps are obtained by using the gas 
particles of the hydro-N-body simulation and are binned in 1'' pixels. The contours correspond to the cluster flux. 
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Figure 2: Spectroscopic temperature map of the simulated cluster of galaxy as derived from the spectroscopic 
analysis of the Chandra "observation" with the package X-MAS 



information carried by gas particles and so usually define physical temperatures. In the years 
a number of temperature functions have been introduced in order to compare hydro-N-body 
simulations with X-ray observations. The most common temperature function used is the so 
called bolometric emission- weighted temperature function defined as 
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where A(T) = CEdE oc VT (see, e.g., pj). To make observations more directly compa- 
rable with simulations we build the software package X-ray Map Simulator (X-MAS) devoted 
to simulate X-ray observations of galaxy clusters obtained from hydro-N-body simulations 
Using X-MAS we clearly show that if the cluster has a complex thermal structure the emission- 
weighted temperature function does not reproduce the spectroscopic estimate and in particular 
it overestimates this valuef^J. This problem is illustrated in Fig. Fig. El Fig. El In fact, in 
the left panel of Fig Q] we show the emission- weighted temperature map of a simulated cluster 
that we use as input for X-MAS. We point to the reader the presence in this map of two shock 
fronts: one in the lower-left corner and the other at the center of the Eastern side of the cluster. 
In Fig. 13 we show the Tspec 

map derived from the data analysis of a 300 ks Chandra ACIS-S3 
X-MAS "observation" of the cluster shown in Fig Q] (see f\ for details). 

In order to make a direct comparison of the "observed" temperature map to the emission- 
weighted one, we decreased the resolution of the latter to match the resolution of the former. 
Thus, in the left panel of Fig. [HI we report the same map of Tgw shown on the left panel of 
Fig. d but re-binned as the map of Fig. [21 To highlight the temperature differences between 
the 

2^spec 9^nd the Tew maps, in the right panel of Fig. El we show (Tew ~ 2^spec)/2^spec- We 
only show the pixels where this difference is significant to at least 3a confidence level, i.e. 
|(Tew — Tspec)/crspec| > 3. This plot clearly shows that there are many regions where the difference 
between Tgw and Tgpec is significant to better than 3a confidence level. Furthermore, for these 






Figure 3: Left panel: emission- weighted temperature map of the simulated cluster of galaxies shown in the left 
panel of Fig. Q] but re-binned to match the spatial resolution of the spectroscopic temperature map shown in 
Fig. 121 Right panel: percentile difference between the spectroscopic and emission- weighted temperature maps. In 
this map we show only the regions where the significance level of the temperature discrepancy is at least 3cr, i.e. 
I (Tew — Tspec)/c"spec| > 3, whcrc (Tspec is the 68 per cent confidence level error associated to Tgpec- 



pixels the discrepancy ranges from 50 per cent to 200 per cent and even more. Of particular 
relevance are two cluster regions showing a shock front in the Tgw map, in the left panel of 
Fig. nj we notice discrepancies of 100-200 per cent, indicating that shock fronts predicted in the 
Tew map are no longer detected in the observed the Tgpec map. 



2 A new formula to estimate spectroscopic temperatures 

If the cluster thermal structure is complex, the Tew may substantially differ from Tgpec- This 
is mainly related to the fact that the projected temperature function is not a good physical 
quantity. From a pure analytic point of view, in fact, the sum of the spectra of two or more 
thermal model is no longer a thermal model Yl- the real world, however, the finite energy 
response and limited energy resolution of the X-ray instruments conspire to distort the observed 
spectra making multi-temperature thermal source spectra fitted by single-temperature thermal 
models which have little to do with the real temperature, but nevertheless are statistically 
indistinguishable from it (see ^ for details). In we show that the temperature that better 
approximates Tgpec is given by what we call the "spectroscopic-like" temperature Tgi. The idea 
behind the derivation of this temperature function is quite simple: if we assume two thermal 
components with densities ni, n2, and temperatures Ti, T2, respectively, requiring matching 
spectra means that 



njCiZ, Ti)^ exp(-^) + njCiZ, Ts)^ exp(- J.; 
^AC(Z,r3i)^exp(-^), 



(2) 



where A is an arbitrary normalization constant and Q{Z^T) is a parametrization function that 
accounts for the total Gaunt factor and partly for the line emission. 

Both Chandra and lOsAM-Newton are most sensitive to the soft region of the X-ray spec- 
trum, so we expand both sides of Eq. I2lin Taylor series, to the first order in E/kT. By equating 
the zero-th and first-order terms in assuming a power-law functional form for T2), and 
extending to a continuum distribution we find that (see for details): 
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It is interesting to note that Tgi weights each thermal component directly by the emission measure 
but, unlike Tew, inversely by their temperature to the power of 3/4. This means that the observed 
spectroscopic temperature is biased toward the coolest regions. To prove the quality of Tgi in 
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Figure 4: As Fig. El but for the spectroscopic-like temperature. 

reproducing Tgpec we do the same test done for Tgw and shown in Fig. El In the left panel of 
Fig. ^ we report the cluster Tq\ map re-binned as the map of Fig. El In the right panel of Fig. 01 
we show (Tgi — Tgpec )/rspec. The presence in this map of fewer pixels clearly indicates that the 
match between Tgi and Tgpec is much better than the one between Tgpec and Tew Furthermore, 
most of these pixels shows very small temperature discrepancies. This demonstrates that the 
Tgi temperature gives a much better estimate of the observed Tgpecthan the widely used Tew 

It is very important to say that, unlike the emission- weighted, the map on the right panel 
of Fig. 01 does not show big discrepancies between Tgi and Tgpec in both shock regions. This 
clearly indicate that Tgi does a much better job than Tew in predicting the projected spectral 
properties of such peculiar thermal features. 

3 Conclusions 

The projected Tgpec of thermally complex clusters is lower than Tew This bias has important 
implications for the study of thermal structures in clusters like shock fronts. In fact in real 
observations shock fronts appear much weaker than what is predicted by emission- weighted 
temperature maps, and may even not be detected. This bias effect is evident on the right panel 
of Fig. n] where we report the spectroscopic-like temperature map of the simulated cluster. We 
notice that the Tgi map appears cooler than the map of Tew Furthermore, both shock fronts, 
which are clearly evident in the Tew map, are no longer detected in the Tgi map. This bias 
may explain why, although numerical simulations predict that shock fronts are a quite common 
feature in clusters of galaxies, to date there are very few observations of objects in which they are 
clearly seen. To conclude we stress that the emission-weighted temperature function may give 
a misleading view of the actual gas temperature structure as obtained from X-ray observations. 
Thus, here we propose to theoreticians to finally discard its use. 
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